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ABSTRACT 


The  purpose  of  this  thesis  was  to  develop  a  technique  by  which  the 
frequency  response  of  a  servomechanism  might  be  rapidly  evaluated  by 
graphical  methods.  The  approach  consisted  of  developing  several  frequency 
response  parameters  related  to  the  location  of  the  poles  and  zeros  of  the 
closed  loop  transfer  function  of  a  servomechanism.  The  systems  considered 
were  the  pure  second  order  (one  pair  of  complex  conjugate  poles) ,  second 
order  with  one  zero  (one  pair  of  complex  conjugate  poles  and  one  real 
zero) ,  and  a  pure  third  order  (one  pair  of  complex  conjugate  poles  and 
one  real  pole).  For  the  latter  two  systems,  the  complex  poles  and  the 
frequency  were  redefined  as  referred  values  relative  to  the  real  pole 
or  zero. 

Loci  of  constant  frequency  response  parameters  were  mapped  onto  a 
complex  plane;  or  loci  of  constant  imaginary  referred  pole  locations  were 
mapped  onto  an  auxiliary  complex  plane,  depending  upon  the  specific 
system  considered.  Using  the  points  obtained  from  the  appropriate 
relationship,  the  frequency  response,  M  and  N  contours,  could  be  rapidly 
and  easily  sketched.  The  M  and  N  contours  could  then  be  used  as  a  basis 
for  rapid  system  evaluation  in  the  frequency  domain. 

The  authors  wish  to  thank  Professor  George  J.  Thaler  for  bringing 
to  their  attention  the  thesis  by  Major  Walter  B.  Patton,  USMC ,  and 
LT  William  B.  Abbott,  III,  USN,  without  which  this  thesis  topic  would 
probably  not  have  been  considered  / 1 / . 

The  authors  are  deeply  indebted  to  Assistant  Professor  Robert  D.  Strum, 
who  acted  as  their  thesis  advisor,  for  the  friendly  support,  encour¬ 
agement,  and  advice  given  so  selflessly. 
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1.  INTRODUCTION 


Each  year  the  importance  and  utilization  of  servomechanisms  and 
automatic  control  systems  increase.  It  naturally  follows  that  more 
papers  are  submitted  pertaining  to  the  design,  analysis,  synthesis, 
and  compensation  of  servomechanisms  and  automatic  control  systems. 

Many  of  these  papers  and  existing  texts  discuss  details  and  laborious 
calculations  which  are  certainly  a  prerequisite  to  analyze  or  synthesize 
the  operation  of  a  servomechanism.  However,  in  many  cases,  a  rapid 
approximation  of  system  operation  would  be  beneficial,  and  possibly 
only  an  approximation  would  be  required  upon  which  an  intelligent 
estimate  could  be  made  on  the  performance  of  the  system  under  study. 

There  are  only  a  few  papers  known  to  the  authors  in  which  the  use  of 
approximations  are  presented  from  which  an  estimate  of  system  performance 
can  be  made,  e.g.,  Floyd's  method  /2/.  One  such  paper  by  Patton  and 
Abbott  /l/  presents  a  method  for  a  rapid  approximation  of  servomechanism 
transient  response. 

A  corollary  should  exist  to  the  Patton-Abbott  development  in  which 
a  rapid  approximation  of  the  frequency  response  of  a  servo  could  be 
determined.  The  use  of  frequency  response  analysis  is  still  an  important 
criterion  in  measuring  servo  performance,  although  the  transient  response 
approach  has  become  very  popular  and  widely  used  in  recent  years,  almost 
to  the  point  of  neglecting  frequency  response.  There  are  many  ways  of 
determining  frequency  response  from  transient  response,  but  lengthy 
calculations  are  necessary. 

It  is  the  purpose  of  this  paper  to  present  a  graphical  means  for 
determining  a  rapid  approximation  of  servomechanism  frequency  response. 
This  approach  eliminates  the  use  of  the  Bode  diagram  and  Nichols' 
chart.  The  M  and  n  contours  will  be  obtained  directly  from  the  closed 
loop  transfer  function  of  the  servo. 
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This  development  will  parallel  the  Pa tton-Abbo tt  approach  so  that 
from  the  same  given  data  (closed  loop  transfer  function,  factored  form) 
it  is  possible  to  obtain  rapidly  the  transient  response  from  their 
paper  and  the  frequency  response  from  this  paper.  This  paper  will 
discuss  in  sufficient  detail  the  development  of  the  frequency  response 
approach  so  that  no  reference  to  the  Patton-Abbo tt  paper  will  be  required. 
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2  GENERAL  DEVELOPMENT 


Since  the  time  domain  or  transient  response  is  used  so  widely  in 
the  initial  phases  of  servo  study,  all  developments  in  this  paper 
originated  with  the  closed  loop  transfer  function  expressed  in  Laplace 
notation  (s-plane)  with  all  initial  conditions  equal  to  zero.  The 
transformed  equation  was  then  converted  into  frequency  response  notation. 
For  this  to  be  a  valid  conversion,  certain  basic  definitions  are  in  order, 
plus  some  additional  restrictions  and  assumptions  used  in  this  devel¬ 
opment.  The  frequency  response  of  a  system  was  defined  in  part  by 
considering  only  steady  state  conditions.  Thus,  all  transient  conditions 
have  vanished.  To  convert  the  transfer  function  from  Laplace  notation 
into  the  frequency  response  equation,  it  was  only  necessary  to  make  a 
substitution  of  the  variable.  This  was  based  on  the  assumption  that 
the  input  was  a  pure  sinusoid  and  that  the  output  was  also  a  pure 
sinusoid  since  the  steady  state  condition  exists.  Furthermore,  this 
forces  the  restriction  that  the  system  must  be  linear.  The  complex 
variable  s  ®  o  +  ju)  then  becomes  the  imaginary  frequency  variable 
since  the  real  part  equals  zero  (a  »  0).  Therefore  s  “  ju);  and  when 
this  variable  substitution  was  made  in  the  closed  loop  transfer 
function,  the  frequency  response  equation  resulted. 

2 . 1  Basic  Definitions  and  Systems  Characteristics 

At  this  point,  a  tabulation  of  the  various  definitions,  restric¬ 
tions,  and  assumptions  would  be  in  order: 

a.  The  system  was  linear. 

b.  The  system  was  stable. 

c.  The  system  was  characterized  by  unity  feedback,  i.e., 

M  =  1  at  oj  =  0. 

do  The  closed  loop  transfer  function  was  known  in  its  factored 

forme 

Based  on  the  above,  some  general  frequency  response  characteristics 
could  be  stated. 
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2 . 2  M-Contour 


a  At  zero  frequency,  M  is  unity. 

b.  At  zero  frequency,  the  slope  of  the  M  contour  is  zero. 

c.  As  the  frequency  approaches  infinity,  M  approaches  zero. 

d.  As  the  frequency  approaches  infinity,  the  slope  of  the  M 
contour  approaches  zero. 

e.  At  the  maximum  or  minimum  value  of  M,  the  slope  of  the  M 
contour  is  zero. 


2 . 3  N-Contour 

a.  At  zero  frequency,  N  is  zero. 

b.  As  the  frequency  approaches  infinity,  N  approaches  ^  tt 

degrees . 

(m  is  the  number  of  zeros  and  n  is  the  number  of  poles  of  the  closed 
loop  transfer  function) 

2.4  Parameters 

To  display  graphically  the  frequency  response,  parameters  were 
selected  to  define  specific  points  on  the  M  and  N  contours.  The 
selected  parameters  were  defined  as; 

a.  M  :  The  value  of  M  where  the  slope  of  the  M  contour  was  zero 

P 

(Does  not  include  the  value  of  M  at  zero  frequency  when  M  ■  1.0). 

b.  00^:  The  value  of  omega  (cu  »  2rrf)  where  the  slope  of  the  M 

contour  was  zero.  (Does  not  include  U)  ■  0) . 


c. 

NP: 

The 

d. 

Mbw: 

The 

e. 

“W 

The 

f. 

NbwJ 

The 

g- 

M,: 

The 

h. 

(U  • 

V 

The 

i. 

N.: 

The 

M^:  The  notation  when  M  ■  1.0,  except  at  ou  ■  0. 
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j.  u>^  :  The  undamped  natural  angular  frequency. 

k.  M  :  The  value  of  M  at  the  undamped  natural  angular  frequency  (U)  ) 

UDn  n 

l.  N  :  The  value  of  N  at  the  undamped  natural  angular  frequency  (<JU  ) 


The  above  selected  parameters  defined  the  general  shape  of  the  M 
and  N  contours.  However,  it  was  discovered  that  with  only  these  points, 
the  M  and  N  contours  were  in  some  cases  too  general;  the  contours  had 
insufficient  shape  to  be  of  significant  value.  Better  approximation  of 
the  M  and  N  contours  were  provided  by  the  selection  of  additional  M,  N, 
and  (JO  points  to  "fill-in11  the  shape  of  the  contours.  It  was  discovered 
that,  mathematically,  the  choice  of  as  the  basis  for  these  additional 
points  reduced  the  M  and  N  equations  to  relatively  simple  relationships 
and  insured  that  the  additional  points  were  located  in  the  critical 
regions  of  the  contours.  The  benefits  gained  by  the  additional  points 
more  than  compensate  for  the  slight  additional  labor  by  providing  a  much 
more  accurate  approximation  to  the  M  and  N  contours.  The  additional 
parameters  selected  were: 

m.  00^:  An  arbitrary  frequency  proportional  to  <u  (u^  02  k  u^)  , 

where  k  was  a  selected  constant  as  discussed  later. 


n. 

V 

The 

va  lue 

of 

M 

at 

V 

o . 

V 

The 

value 

of 

N 

at 

A  sketch  of  the  general  M  and  N  contours  is  shown  in  Fig.  2-1. 


2 . 5  Classes  of  Servomechanisms  Studied 

From  Fig.  2-1,  it  can  be  seen  that  it  was  relatively  easy  to 
construct  an  approximate  frequency  response  if  the  above  defined  para¬ 
meters  were  known.  This  paper  contains  the  above  parameters  presented  in 
graphical  form.  This  paper  also  contains  the  technique  utilized  in  the 
parameter  determination  such  that  the  frequency  response  can  be  rapidly 
approximated  for  specific  classes  of  systems*  The  three  systems  covered 
in  this  paper  were: 

I.  Pure  Second  Order  -  The  closed  loop  transfer  function  consisted 
of  two  complex  poles  only. 
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Fig.  2-1  General  M  and  N  Contours  Showing  Parameters  Chosen 
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II  Second  Order  with  One  Zero  -  The  closed  loop  transfer  function 


consisted  of  one  pair  of  complex  poles  and  one  real  zero. 

III.  Pure  Third  Order  -  The  closed  loop  transfer  function  consisted 
of  one  pair  of  complex  poles  and  one  real  pole. 

2  6  Presentation  of  Parameters 

The  graphical  presentations  for  the  pure  second  order  system  were, 
for  the  most  part,  loci  of  constant  values  of  the  various  parameters 
plotted  on  the  s-plane.  For  the  other  two  systems  covered  in  this 
paper,  it  was  necessary  to  express  the  coordinates  of  the  complex  poles 
as  a  function  of  the  magnitude  of  the  real  pole  or  zero.  In  accomplishing 
this,  the  angular  frequency,  U),  also  became  a  function  of  the  real  pole 
or  zero  location  Thus  the  magnitude  and  phase  shift  for  the  frequency 
response  was  shown  to  be  a  function  of  the  complex  pole  location  relative 
to  the  real  pole  or  zero.  The  graphical  presentations  for  the  second 
order  plus  one  zero  and  for  the  pure  third  order  systems  were  plotted  in 
two  ways.  First,  by  loci  of  constant  values  of  the  various  parameters 
plotted  on  a  referred  complex  plane,  that  is,  with  ordinate  and  abscissa 
being  the  complex  pole  location  relative  to  the  real  pole  or  zero. 

Second,  by  loci  of  constant  referred  imaginary  part  of  the  complex  pole 
plotted  against  referred  real  part  of  the  complex  pole  as  abscissa  versus 
the  chosen  parameter  as  ordinate. 

2 . 7  Digital  Computer  Utilization 

The  solutions  of  the  equations  that  will  be  presented  later  when 
a  specific  system  is  discussed  were  obtained  from  a  Control  Data  Corp¬ 
oration  high-speed  digital  computer  Model  1604,  located  at  the  U.  S.  Naval 
Postgraduate  School,  Monterey,  California  (Appendix  A).  Without  the 
availability  of  a  high  speed  digital  computer,  it  is  doubtful  that 
this  paper  could  have  been  developed  without  the  utilization  of  an 
analog  computer  plus  excessive  labor. 
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3  THE  PURE  SECOND  ORDER  SYSTEM 


The  pure  second  order  system  treated  in  this  section  was  pre¬ 
viously  defined  as  a  system  whose  closed  loop  transfer  function  consisted 
of  only  two  complex  conjugate  poles.  This  system  has  been  extensively 
treated  in  most  servo  textbooks,  usually  in  terms  of  the  damping  ratio 
and  the  undamped  natural  frequency.  Here,  however,  it  will  be  treated 
by  the  rectangular  location  of  the  closed  loop  poles.  The  closed  loop 
transfer  function  of  the  pure  second  order  system  was  stated  as  /!/: 


CL  j-  b 

(SiC’f-t  <7 


(3.1) 


The  poles  of  this  system  were  plotted  on  the  s-plane  in  Fig.  3-1. 
Converting  Eq.  (3.1)  to  the  frequency  response  form,  it  became 


3 . 1  M  and  (JO  Calculations 

Equation  (32)  rearranged  became 


from  which 


(3.2) 


(3.3) 


M 


CL 


i 

fo 


I 


UL* 


(3.4) 


Note  that  if  the  system  gain  were  constant,  then  "a"  and  nbn  would  be 
fixed  and  the  magnitude  would  be  a  function  only  of  system  frequency, 
i .  e  .  ,  M  =  f  ( uo) 
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Fig.  3*1  Poles  of  the  Pure  Second  Order  System  (s-plane) 
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3  2  M  and  cd  Determination 

£2 

CD  was  obtained  by  finding  the  value  of  CD  where  the  slope  of  the 
P 

M  versus  CD  curve  was  equal  to  zero.  For  mathematical  convenience,  this 
was  done  by  first  squaring  Eq .  (3.4),  then  taking  the  first  derivative 
and  setting  the  numerator  equal  to  zero. 


(a  -<  b  )  2.  (d  -f  b  —u/')(-2.co)  +  8  o.'co  J  —  0 


(3.5) 


Simplifying,  this  becomes 

Us  ~  Usp  —  ~\j  jo  —  CL 
and  CD  =  0  (see  Section  2,4b) 


(3.6a) 

(3.6b) 


Substituting  Eq .  (3.6a)  in  Eq . 


(3.4)  gave  the  value  of  M  at  (D 


*  I** 

2.ob 


(3.7) 


Note  that  from  Eq  (3.6),  M  would  not  occur  whenever  "a"  was  less  than 

P 

Mbn,  i.e.,  when  the  real  part  of  the  pole  was  greater  than  the  imaginary 
part  (Figs,  3-2  and  3-3) 

3.3  CD,  Determination 
bw 

The  equation  for  the  bandwidth  frequency  was  obtained  by  substituting 
the  bandwidth  magnitude  (M  =  .707)  into  Eq .  (3.4)  and  solving  for  the 

frequency  This  resulted  in 

*■  y/'*~v+ um'/ai4  <3-s> 


(Fig  3-4) 
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3.4  M  and  uu  Determination 
cun  n 


The  undamped  natural  frequency  was  defined  as 


^  y  a.  b 

This  equation  expressed  u)  directly  in  terms  of  "a"  and  "b". 
The  substitution  of  into  Eq .  (3.4)  yielded 


(3.9) 
(Fig.  3-3) 


(3.10) 


Equation  ( 3. 10)expressed  M  in  terms  of  the  closed  loop  poles  of  the 

ctn 

system.  (Fig  3-6) 

3  5  Determination 

As  defined  previously,  0)^  was  the  frequency,  other  than  zero,  at 
which  M  was  unit?/.  Letting  M  =  1*0  in  Eq .  (3.4)  and  solving  for 
yielded 


(jL\  — 


(3.11) 


Simplifying  this  expression  by  recalling  zzr. 1 

resulted  in  OJf  ^  ^f'Z*  CAJy,  (3.12) 

(Fig .  3-7) 

3  6M,  and  uo  Determination 
 k_ k 

As  stated  in  Section  2,  intermediate  values  of  M  and  N,  in  addition 
to  those  determined  above,  would  result  in  a  smoother  plot  approximating 
the  M  and  N  contours.  The  values  of  k  chosen  for  this  class  of  servo 
were  .5  and  1.5  These  particular  values  were  selected  as  being  those 
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points  which  would  be  most  beneficial  in  smoothing  the  plot  of  the 
M  and  N  contours,  and  because  the  equation  was  the  same  for  k  =  .5 


as  tor  k  =  1 «  5 

In  Section  2, 

CO L  CC  k  LOy 

or 

H, "  {K  =  i~5  UJ? 

and 

~  ^°p  —  L5 

(Fig  3-7) 

The  value  of  U) 

and  U)  were  substituted  into 

ki 

k2 

2 


1. 2.  CLAh' 


(.%  '4 


f  i/f 


(3.13) 


Where  ,  recalling  CLy} ^  y  cC  „  Loci  of  constant  values 

of  were  plotted  in  relation  to  "a"  and  nbn.  (Fig.  3-8) 


3. 7  N  Calculations 
From  Eq .  (3-3) 


r /.■  a 

K  L  \oj)  — 

^  o 


...  -t  h  -co  c  'LO  lOj 


+Zalc 


—  £  M  ^ 


:N 


from  which, 


A/ ~  -  oaclan 


to  oj 


'x  lA  2, 

cl  - tb  -  J 


(3.14) 


Equation  (3.14)  was  then  the  basic  N  equation  for  the  pure  second  order 
system. 
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3.8  N  Determination 

J2 

u>p  was  previously  calculated  by  Eq .  (3.6).  The  substitution  of 
Eq .  (3.6)  into  Eq,  (3.14)  resulted  in  the  solution  of  related 
directly  to  "a"  and  "b",  the  closed  loop  poles  of  the  system. 


—  c.  ret  an 


(3.15) 


Loci  of  constant  values  of  N  were  then  determined  for  all  possible 

P 

"a"  and  "b"  values.  (Fig.  3-9) 


It  is  pertinent  to  note  that  the  equations  presented  were  in  most 
cases  not  solved  as  stated.  Algebraic  manipulations  were  performed  to 
make  the  equation  compatible  to  digital  calculations  by  writing  the 
equations  in  terms  of  "b"  as  a  function  of  "a"  and  the  particular 
pa rameter . 


3.9  N  Determination 
(On 


N  was  determined  in  a  manner  similar  to  N  .  The  value  of  <0  in 
(On  p  n 

terms  of  "a”  and  nbM  was  stated  in  Eq .  (3*10).  This  relationship  of 
to^  was  substituted  in  Eq,  (3.14)  and  resulted  in 


[\J  =  —  iU'cia  n 

Vm 


1 


oAcl  ■ 


L4  b 


L  &- 


a.  4 


i-  An 
h  -ha 


- 


which  simplified  to 


-arc-tan 


f  laXa-t  hl  l 

L  0  \ 


(3.16) 


or 


N 

ton 


“90 


o 


Thus,  for  a  pure  second  order  system  was  a  constant  and  equal  to 


-90 


-13- 


3  10  N,  Determination 
bw 

The  equation  for  in  terms  of  "a"  and  Mbn  was  obtained  by 

substituting  Eq .  (3.8)  into  Eq .  (3.14).  This  resulted  in 


N,  —  -arcian 

b w 


f ZcC'ltf- c?+ 1  AM ict+  b<L  '  1 

L  Zcl-  IA\4  i/o4 4  b4  ’ 


(3.17) 


(Fig.  3-10) 


3.11  Determination 

Equation  (3.12)  determined  the  values  of  U)  when  M  =  1.0  for  non-zero 
frequencies.  When  the  expression  of  (D^ ,  determined  from  Eq .  (3.12), 

was  substituted  into  Eq,  (3.14), 


li 


<xrcl 


dot 


n 


2,  CL 


a 


ocC 


-  bx 


(3.18) 


Solutions  to  Eq.  (3.18)  were  plotted  as  constant  loci  for  all  values 
of  nan  and  Mb".  (Fig.  3-11) 


3.12  N.  Determination 
k 


As  stated  above,  the  k  values  selected  were  .5  and  1.5  where 
U)^  ~znJk  OJ^p  #  In  a  manner  similar  to  the  determination  of  N  above, 


(3.19) 


from  which  values  were  determined  explicitly  in  terms  of  “a"  and  nbn 
for  both  values  of  k.  (Figs.  3-12  and  3-13) 


3 .  1 3  Use  of  the  Curves 

At  this  point,  all  the  selected  parameters  necessary  to  sketch  the 
M  and  N  contours  have  been  determined  from  the  location  of  the  closed 
loop  poles.  With  the  values  of  the  closed  loop  poles  known  and  the 
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curves  plotted  as  described  above,  a  rapid  graphical  approximation  of 

the  frequency  response  of  the  pure  second  order  system  can  be  constructed. 

For  convenience,  the  following  recommended  procedure  is  presented  for 

the  use  of  the  pure  second  order  curves: 

a.  Enter  Fig.  3-2  with  the  complex  pole  coordinates  -a  +  jb,  and 

obtain  M  *  Enter  Fig.  3-3  with  -a  +  jb,  and  obtain  U)  .  (Note  that  U) 

P  P  P 

is  the  ordinate  in  Fig.  3-3.) 

b  Similarly  obtain  (1)  (M  =  .707)  from  Fig.  3-4,  0)  and  M 

bw  n  ODn 

from  Figs,  3-5  and  3-6  respectively. 

c.  Obtain  (M  =  1.0)  from  Fig.  3-7,  entering  with  the 
determined  above* 


wi 


d.  Obtain  for  both  =  .  5  and  =  1.5  from  Fig.  3-8,  entering 
th  "a"  and  "b".  Determine  U)  ,  (a CO^)9 


and 


uo,  /  m.  - 

k2  *  v 


'\li. 


5  cop  , 


from  Fig.  3-7,  entering  with  the  value  of  obtained  from  step  (a)  above. 

The  M  contour  for  the  specified  closed  loop  poles  may  then  be 

sketched,  having  easily  obtained  the  M  and  corresponding  U)  for  six  known 

points  and  realizing  that  M  s  0  and  has  a  slope  of  zero  at  (D  =  0,  and 

also  that  the  slope  is  zero  at  (JO  and  at  (D  =  °°. 

P 

The  N  contour  points  are  obtained  in  the  same  manner  by  using 

Figs.  3-9  through  3-12  to  obtain  the  values  of  N  corresponding  to  the 

same  frequencies  found  for  the  M  contour  above.  The  one  exception  is 

N  ,  which  is  always  equal  to  -90°  for  this  system.  Using  these  N, 
ci)n 

0)  values  and  knowing  that  N  =  0  and  has  zero  slope  for  0)  =  0 ,  and  also 
that  N  =  -180°  for  (i)  =  00  for  a  pure  second  order  system,  the  N  contour 
can  now  be  accurately  sketched. 


Note  that  this  paper  does  not  present  curves  which  will  give  any 

values  at  frequencies  higher  than  U)  ,  since  it  is  unlikely  that  there 

bw 

would  ever  be  any  interest  in  the  region  beyond  the  bandwidth  frequency. 
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3  1 4  Exarrple 

Given:  -a  +  jb  =  -1  +  j2 

Solution: 

SPECIFIC  DATA 

TABLE  OF  PARAMETER  VALUES 


X 

(JO 

(JO  Fig.  No. 

M 

M  Fig.  No. 

N 

N  Fig.  No. 

X 

x 

X 

X 

X 

X 

p 

If! 

3-2 

1.Z5 

3-3 

3-9 

bw 

2. 9* 

3-4 

.  707 

N.A. 

3-10 

(JO 

n 

z.zS' 

3-5 

1. 17 

3-6 

-90° 

N.A. 

1 

Z-4-Z 

3-7 

1.0 

N.A. 

3-11 

ki 

17  2 

3-7 

/  i/~r 

3-8 

-o^r° 

3-12 

k2 

x.io 

3-7 

lXjf 

3-8 

-82° 

3-13 

GENERAL  DATA  (for  all  pure  second  order  systems) 


(JO 

M 

dM/du) 

N 

dN/d(JU 

0 

1.0 

0 

0° 

0 

(JO 

M 

0 

P 

P 

00 

0 

0 

-180° 

0 

Fig  3-14  (page  17  )  compares  the  frequency  response  obtained  from 
this  rapid  approximation  method  to  the  actual  M  and  N  contours  obtained 
from  a  digital  computer  program  of  the  general  M  and  N  equations. 

Note :  Although  a  scale  is  shown  on  these  pure  second  order  curves  for 

convenience,  the  loci  are  really  independent  of  scale  as  they  depend 
only  on  the  ratio  of  "a"  to  1 1  b,f  Hence,  any  convenient  scale  for  o9 

and  (JO^  may  be  substituted  for  the  scales  shown  as  long  as  all 
scales  are  in  agreement 
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n:rj: 


By  Approximating  Sketch 

fH" 


Calculated  Curve 

SfiilSMSH8*S 


Fig.  3-1^ 

Example (Frequency  Response 
for  Pure  Second  Order  System) 
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L8 


19 


■21- 


-22- 


-2  3- 


“25 


B 


|::i 


Fig.  3-10 

N^w  for  Pure  Second  Order  System 
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4.  THE  SECOND  ORDER  SYSTEM  WITH  ONE  ZERO 


The  addition  of  a  real  zero  to  a  pure  second  order  system 
resulted  in  a  system  whose  closed  loop  transfer  function  was 


% 

i 


■;i  4  If  (s  f-d) 


(4.1) 


The  singularities  of  this  system  were  plotted  on  the  s~plane  in 

2  2 

Fig  4-1.  The  constant  a  +  b  ,  usually  interpreted  as  system 

d 

gain,  reflected  the  assumption  of  unity  feedback  as  stated  previously. 

By  letting  s  =  jd),  the  frequency  response  form  of  Eq.  (4.1) 
resulted  in 


(4.2) 


This  expanded  and  rearranged  became 

(chi  i  d^j -4  j  (gjCL  +  lo 

(cla-i  cl  .o  -  cl  to)  +  j  (z  daw) 


(4.3) 


Dividing  the  numerator  and  denominator  of  Eq 


(4.3)  by  d3, 


(4.4) 


From  Eq.  (4.4),  it  was  seen  tha t  "a",  "b",  and  uo  may  all  be  referred  to  the 
zero  location,  d. 
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s  -  plane 


Fig.  4*1  Poles  and  Zero  of  the  Second  Order  System  with  One  Zero 
(s- plane) 
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From  Eq ,  (4.4), 


(4.5) 


(4.6) 


and 


r\  £  40 

x  M| 


(4.7) 


These  quantities  were  illustrated  in  Fig.  4-2.  Equation  (4.4)  then 
became 


Thus  the  frequency  response  for  the  second  order  system  with  one  zero 
became  a  function  of  the  "referred"  location  of  the  complex  roots  and 
the  "referred"  frequency.  This  means  then  that  the  frequency  response 
was  the  same  for  all  systems  having  the  same  relative  pole-zero  con¬ 
figuration  when  the  frequency  was  also  referred  to  the  zero  location. 
(Fig.  4-3) 

4. 1  M  and  Cl  Calculations 

Equation  (4.8)  was  expressed  in  terms  of  its  magnitude  and  phase 
angle  by 
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jdj 


m-plane 


o 

-1 


m. 


■>•  y>ij 


X 


Fig  4-2  Referred  Quantities  of  and  n>^,  Second  Order  System 
with  One  Zero 
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System  A- 


System  B: 


a  =  1.0 


a  =  ,01 


b  =  2.0 
d  =  1.0 


b  =  .02 
d  =  .01 


ml  /d/~ 


1.0 


m, 


/d/ 


1.0 


m. 


‘  /d/  2< 


m 


2 


b 

=5  — 

/d/ 


2.0 


0 

0 


1 


01 


2  3 

.02  .03 


4  <oA  (d  =  1.0) 

.04  (i)  (d  =  .01) 
B 


Fig.  4~3  Frequency  Response  of  Two  Systems  with  the  Same  Relative 

Pole-Zero  Configuration,  Second  Order  System  with  One  Zero 
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From  this, 


V  (ml  4-  si  (inj+yrty' 

"V 6?lf -f m*  -ii1 )  +~4infj^ 


(4.10) 


Equation  (4.10)  was  the  basic  M  equation  for  the  second  order  system 
with  one  zero. 


4.2  M  and  Cl  Determination 

_JL _ £ _ 

The  equation  for  Cl ^  was  obtained  by  first  squaring  Eq.  (4.10), 
then  differentiating  with  respect  to  Ci  and  setting  the  numerator  of 
this  equal  to  zero.  This  reduced  to 


SLf  =  +7  1  (m*  -  mf)  + 1 


(4.11) 


Substituting  this  into  Eq.  (4.10)  and  simplifying, 


MP= 


1 


7nj±Z  7 nf  wfc  +  yn£ 


>c"y  ■mj  +lm,mz  +  mz~Z7n,+Z?nl  +  l  -hZ(ml-)>ix-l) 


(4.12) 


Loci  of  constant  values  of  m2  were  solved  by  digital  computer  and  were 
plotted  against  m,  versus  M  and  Cl  .  (Figs.  4-4  and  4-5). 

1  p  p 


4.3  0,  Determination 

DW 

The  equation  for  the  referred  bandwith  frequency  was  obtained  by 
substituting  M  *  .707  into  Eq.  (4,10)  and  solving  for  Q^.  After 
considerable  algebra,  this  reduced  to 


(Fig.  4-6) 


M  f  +  mx  + 


[l(rn.t+mC)%  [ 


7n( 


(4.13) 
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4.4  MU  and  Q  Determination 
Un  n 

The  referred  undamped  natural  frequency  was  defined  for  this 
system  as 


<L 


(4.14) 


This  expressed  Cl  as  a  function  of  m,  and  rru .  (Fig*  4-7) 
n  1  l 

Substituting  Eq .  (4.14)  into  Eq  (4*10)  and  simplifying  yielded 


(Fig „  4-8) 


(4.15) 


4*5  Cl ^  Determination 

As  defined  previously,  Cl ^  was  the  referred  frequency  other  than 
zero  at  which  M  is  unity*  The  equation  for  Cl^  was  found  by  letting 
M  38  1*0  in  Eq  (4.10)  and  solving  for  . 


-ft- 1  -  V  (ia 


(4.16) 


(Fig.  4-9) 

4.6  and 

The  use  of  the  curves  mentioned  up  to  this  point  for  this  class 

of  system  provided  the  four  main  points  required  to  sketch  roughly  the 

M  contour.  In  order  to  "fill -in'*  the  blank  portions  of  the  contour  on 

both  sides  of  fi  ,  it  was  decided  to  provide  curves  for  three  additional 
P 

points  on  the  M  contour  in  the  same  manner  as  was  done  for  the  pure 
second  order  system.  The  use  of  these  additional  points  was  not 
necessary  to  sketch  the  frequency  response;  however,  their  use  provided 
a  more  accurate  approxima t ion  of  the  M  contour 


°k 


Determina  tion 
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As  in  the  pure  second  order  system,  the  following  definition  was  made 


where  k  is  a  constant 


Substituting  into  Eq  (4  11)  for  Q  , 


(4.17) 


SLk  =  T/  -k  i  k f  (jn^ls  -  Z  (rf-yni)  +  1 


(4.18) 


Values  of  .5,  1.5  and  2,5  were  selected  for  k.  Mt  solutions  were 

k 

obtained  for  each  k  value  by  solving  Eq.  (4.18)  for  and  then  sub¬ 
stituting  this  value  into  the  basic  M  equation  with  the  same  values  of 
m^  and  m^  This  was  accomplished  handily  by  one  digital  program. 

(Figs  4-10,  4-11 ,  4-12,  and  4-13) 

4.  7  N  Calculations 

Referring  to  the  basic  closed  loop  transfer  function  in  frequency 
form,  Eq.  (4.9),  the  general  phase  angle  equation  for  the  second  order 
system  with  one  zero,  was  written  as 


N=  ar?l 


Si  (jW  j  Ylj  4“  111  ^  — XL  ) 


■O  Zrn.SX  V  m  S -St 


(4.19) 


The  previously  obtained  parameter  frequency  equations  (Q  ,  £1  ,  Q  , 

p  bw  n 

,  and  the  three  Q^'s)  were  solved  by  digital  computer  programs.  The 
resulting  parameter  frequencies  were  then  substituted  into  the  general 
N  equation  (4.19)  along  with  the  corresponding  values  of  m^  and  m^ .  The 
following  figures  resulted,  all  plotted  with  constant  m^  loci  against  m^ 
versus  N: 


Fig  4-14  N  as  function  of  m-  and  nu 
d  12 


Fig 

4-15 

Nbw 

as 

function  of 

m^  and 

m2 

Fig 

4-16 

N0n 

as 

function  of 

and 

m2 

Fig. 

4-17 

as 

function  of  m ^  and  m^ 

Fig 

4  18 

ni 

ki 

as 

function  of 

m^  and 

m^ 

for 

kl 

=  .5 

Fig 

4-19 

N, 

k2 

as 

function  of 

m^  and 

for 

k2 

=  1.5 

Fig. 

4-20 

N, 

k3 

as 

function  of 

m  ^  and 

m2 

for 

k3 

=  2.5 
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4 . 8  Example 

Given:  Complex  poles;  ~a  +  jb  =  -2.0  +  j  2.0 

Real  zero;  -d  =~2.0 


Solution: 

SPECIFIC  DATA 

TABLE  OF  PARAMETER  VALUES 


GENERAL  DATA  (for  all  Second  Order  plus  One  Zero  Systems) 


Cl) 

M 

dM/du) 

N 

dN/duu 

0 

10 

0 

0° 

tt) 

M 

0 

-- 

p 

P 

00 

0 

0 

p 

VO 

o 

o 

0 

Fig  4-21  (page  39  )  compares  the  frequency  response  obtained  from 

this  rapid  approximation  method  to  the  actual  M  and  N  contours 
obtained  from  a  digital  computer  solution  of  the  general  M  and  N 
equations » 
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bw  for  Second  Order  System  with  One  Zero 


| 
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5  THE  PURE  THIRD  ORDER  SYSTEM 


The  addition  of  a  real  root  to  the  pure  second  order  system 
resulted  in  a  new  system  which  will  be  referred  to  as  the  pure  third 
order  system  The  closed  loop  transfer  function  of  the  pure  third 
order  system  had  one  pair  of  complex  conjugate  poles  and  one  real  pole. 

The  closed  loop  transfer  function  for  the  unity  feedback  third 
order  system  was  written  as 


(5.1) 


Fig  5-1  defined  these  roots  on  the  s-plane. 


Equation  (5.1)  was  simplified  as  follows: 


(5  2) 


and  expanded  became 


tl  4*j  4-  Z  tie  4  If)  4  cf nf+If) 


(5.3) 


Letting  s  =  ju)  in  Eq  (5,3)  resulted  in  the  frequency  response  form 
of  the  transfer  function 


uf('X -i-Zac-  *■)*) •}  C  (of  •- Hf)  (5 •  4) 


or 
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Fig  5-1  Poles  of  the  Pure  Third  Order  System  (s-plane) 
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Numerator  and  denominator  of  Eq  (5.5)  divided  by  c  : 


(5.6) 


11  ll 

Observing  that  a  power  of  c  exists  in  all  terms  corresponding  to 
the  same  power  of  "a",  "bn  and  (A);  the  following  definitions  were 
made : 


(5.7) 


rh 


A 


il 


A  co 


(5.8) 

(5.9) 


These  quantities  were  illustrated  in  Fig.  5-2.  Equation  (5.6)  then 
became 


-f'  ?C-y 


% 


jf/ij  -i  il'x~SL  ("Xil,^  ‘ij  4  i  '-JX  J 


A 

=  Me- 


(5.10) 


n^  and  n^  were  the  complex  pole  locations  referred  to  the  real  pole 
location  and  G  was  the  referred  frequency.  The  frequency  response 
for  the  pure  third  order  system  was  then  a  function  of  the  referred 
complex  pole  and  the  referred  frequency.  Thus  all  systems  having  the 
same  relative  pole  configuration  will  have  identical  response  if  the 
frequency  were  also  referred  to  the  real  pole  location. 


5 „ 1  M  and  G  Calculations 

The  magnitude  portion  of  Eq .  (5.10),  after  algebraic  manipulation, 
reduced  to 


(5.11) 
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Fig.  5-2  Referred  Quantities  and  n^ ,  Pure  Third  Order  System 
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Equation  (5.11)  was  the  basic  M  equation  for  the  pure  third  order 
system. 


5 . 2  Discussion  of  the  M  Contours 

It  was  found  that  the  M  contour  for  a  pure  third  order  system 

will  be  one  of  three  basic  "shapes"  as  shown  in  Fig.  5-3a,  b  and  c, 

depending  on  the  values  of  n^  and  n^*  Figure  5- 3a  was  the  situation 

when  there  was  no  peak.  No  value  of  G  was  obtained  and  hence  no  values 

P 

for  G^  or  the  G^'s.  Figure  5-3b  was  similar  to  the  second  order  M 
contour  shape  and  had  a  single  peak.  For  this  general  "shape"  points 
were  obtained  for  G^,  G^,  G^,  G^ ,  G^  SG^  ,  and  G^  . 


The  third  general  "shape"  is  Fig.  5-3c,  which  resulted  from  the 

G^  equation  (5.12)  which  yielded  two  non-zero  frequencies  at  which 

the  M  contour  slope  was  zero,  one  being  a  "peak"  frequency  and  the 

other  being  a  "minimum"  frequency.  This  latter  frequency  was  designated 

G  to  avoid  being  confused  with  the  G  .  It  was  also  discovered  that 
m  p 

with  this  shape  of  M  contour,  it  was  possible  to  have  M  less  than 

P 

unity.  Another  possibility  was  to  have  less  than  M^w  and  at  the  same 
time  >1^  greater  than  one.  This  resulted  in  two  values  of  G^,  two  values 
of  Gk  ,  Mk  ,  Gk  ,  and  (one  relative  to  and  the  other  relative  to 

M  ),  and  three  values  of  G,  (one  for  each  crossing  of  the  M  =  .707  line), 
m  Dw 

in  addition  to  the  usual  points  at  G  ,  G  ,  and  G  . 

p  m  n 

Figure  5-4  might  be  thought  of  as  the  complex  n-plane  and  shows  the 
regions  in  which  certain  n^  and  n^  combinations  produce  one  of  the 
three  above  M  contours.  This  figure  indicates  the  general  shape  of  the 
M  contours  that  result  from  the  known  n^  and  n^  values. 


5.3  M  and  G  Determination 
— E_ E 

The  equation  for  G  was  obtained  by  squaring  Eq .  (5.11),  then 
differentiating  with  respect  to  G  and  setting  the  numerator  of  this 
equal  to  zero.  This  reduced  to 
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Fig  5-3  The  Three  Basic  M  Curve  Shapes  for  Pure  Third  Order  System 


(a) 


I 


The  plus  and  minus  values  of  the  inner  square  root  quantity  were  both 
valid  and  resulted  in  two  possible  solutions  for  G^.  The  consequences 
of  two  solutions  were  discussed  in  detail  in  paragraph  5.2  above. 

Loci  of  constant  values  of  n^  were  solved  by  a  digital  computer 

program  and  were  plotted  against  n_  versus  G  and  G  in  Figs.  5-5a 

1  pm 

and  5-5b  In  the  same  computer  program,  values  of  were  obtained 

for  corresponding  values  of  n^,  G^,  and  n^  by  substituting  these 

known  quantities  into  the  basic  M  equation,  Eq .  (5.11).  Figures  5-6a, 

5~6b,  and  5~6c  contain  M  plotted  as  a  function  of  n.  and  n0 .  No 

P  1  ^ 

specific  M  equation  was  derived  due  to  the  obvious  complexity 
involved  in  substituting  Eq.  (5.12)  into  Eq.  (5.11)  for  G. 


5,4  G,  Determination 

DW 

The  equation  for  the  referred  bandwidth  frequency  was  obtained 
by  substituting  M  =  .  707  into  Eq  (5.11)  and  solving  for  n^  as  a  function 

2 


of  Gl  and  n 

DW 


This  involved  extensive  algebra  which  reduced  to 


72  —  L  4*  r ^ w /j. 4  a  ?  \  ~ 

1  U4rf”i  j~[l  1  J 


■±, 


(5.13) 


The  +  signs  were  both  valid  here  also  and  their  effect  will  be  discussed 
later.  (Fig.  5-7) 

5.5  MU  and  G  Determination 
k2n  n 

The  undamped  natural  frequency  for  the  pure  third  order  system  was 
defined  as 


n?+7Lx 


^8 


*«  * 

G^  was  expressed  as  a  function  of  n^  and  n9 


(5.14) 


(Fig.  5-8) 


Substituting  this  into  Eq.  (5.11)  and  simplifying  yielded 


>  -*  1 7:-T 


V  4 n. Oh  +*£  XK -f  l) 


(5.15) 


was  thus  expressed  as  a  function  of  n^  and  n£-  (Figs.  5-9a  and  5-9b) 
5.6  0^  Determination 

The  referred  frequency,  other  than  zero,  at  which  M  ■  1.0  was 
defined  as  Q  . 

Letting  M  *  1.0  in  Eq.  (5.11)  and  solving  for  resulted  in 


J), 


n}4 


r- 


±1% 


n, 


-4n?nZ+. 


2 .5 


(5.16) 


(Figs.  5-10  a,  b,  c,  and  d) 

5.7  and  0^  Determination 

For  the  same  reasons  as  discussed  in  Section  2.4,  it  was  found 
desirable  to  obtain  "fill-in"  points  on  both  sides  of  fi  .  This 
was  done  by  defining 

=  kQp  (5.17) 

which  eliminated  the  square  root  of  k  as  used  for  the  previous  two 
classes  of  systems. 

The  values  of  k  chosen  were  .8,  1.4,  and  1.8;  therefore 

Q,  *  8  Op 

kl 


Cl  ■  1.4  Qp 
k2 


a  ■  i  .  8  r*p 
k3 


(5.18) 

(5.19) 

(5.20) 
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Fig.  5-11  was  the  almost  trivial  plot  of  these  three  equations, 
included  for  optional  use  in  determining  the  0^  values  once  has 

been  determined. 

Loci  for  the  three  M^’s  were  calculated  by  a  digital  program  which 
utilized  Eqs.  (5.11),  (5.12),  (5.18),  (5.19),  and  (5.20);  thus  it  was 
not  necessary  to  attempt  to  obtain  any  specific  equations. 

(Figs «  5* 12  a , b ,  5-13  a,b,  5-14) 


5 . 8  Determination  of  N  Curves 

Referring  to  the  basic  closed  loop  transfer  function  in  frequency 
form,  Eq .  (5.10),  the  general  phase  angle  equation  for  the  pure  third 
order  system  was  written  as 


/v  =  -  aKian f— tznill.  L 

L  ~~a  (%  1 )  -I  7  if* -j  71% 


(5.21) 


The  previously  obtained  parameter  frequency  equations  (fi^, 

,  and  Q^)  were  solved  by  digital  computer  programs.  The  resulting 
parameter  frequencies  were  then  substituted  into  the  general  N  equation, 
Eq,  (5.21),  along  with  the  corresponding  values  of  n^  and  n^ •  The 
following  figures  resulted,  all  plotted  with  constant  n^  loci  against 
n^  versus  N: 


Fig .  5-15  a 
Fig .  5-15  b 
Fig,  5-16  a 
Fig.  5-16  b 
Fig  5-17  a 
Fig  5-17  b 
Fig  5-17  c 
Fig,  5-18  a 


N  as  a  function  of  nn  and  n« 

P  1  2 

N  as  a  function  of  n.  and  n0 

m  12 

N.  as  a  function  of  n,  and  n0 
bw  12 

^bw  (extenc*ec*)  as  a  function  of  n^  and  n 
as  a  function  of  n^  and  n^ 

(expanded)  as  a  function  of  n^  and  n^ 
(2nd  value)as  a  function  of  n^  and  n^ 
for  maximum  (fi^) 
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Fig. 

5-18 

b 

\ 

for 

minimum 

<nm> 

Fig. 

5-19 

a 

Nk 

k2 

for 

ma  x  imum 

<V 

Fig. 

5-19 

b 

Nk 

*2 

for 

minimum 

<fi„> 

Fig. 

5-20 

Nk 

for 

maximum 

<v 

An  equation  for  was  obtained  by  substituting  Eq .  (5.14)  into 
the  basic  N  equation,  Eq .  (5.21), 


fi-n  =  -  anian 

—  —  o.rctan  >  -7  J - -  I 

n  *  -J 


(5. 


Loci  of  were  plotted  as  a  function  of  n^  and  r\2  •  (Fig.  5-21) 


22) 
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5 . 9  Example 


Given  Complex  poles*  -a  +  jb  =  -.8  +  j  4.0 


Real  pole: 


-c  =-2.0 


Solution  : 
SPECIFIC  DATA 


nl  =  /c/  =  -4 

n2  =  777  =  2-° 


TABLE  OF  PARAMETER  VALUES 


X 

G 

X 

(JU  kirk 

X 

G  Fig. No. 

X 

M 

X 

Mx  Fig. No. 

N 

X 

N  Fig. No. 

X 

p 

/  8S 

5. 7© 

5- 5a 

//f 

5- 6a  ,b 

-  10c c 

5- 15a 

m 

l.t% 

5-  5b 

7  6/ 

5- 6c 

06,0 

5-15b 

bw 

7. 

*  6© 

5-7 

-  707 

N.A. 

-172- O'1 

5-  16b 

bw 

N  A 

— — 

— 

.  707 

N.  A. 

N.A. 

— 

bw 

N .  A . 

— 

.  707 

N.A. 

N.A. 

n 

2.0* 

4-CS 

5-8 

U2. 

5-  9b 

-JO/ 

5-21 

1 

Z.14- 

7,7.6 

5-10a  ,b 

1.0 

N.A. 

-l6£  t° 

5-17a  ,b 

1 

I.zz 

72/0 

5- 10c ,  d 

1.0 

N.A. 

-£S° 

5-  17c 

ki  * 

1,02 

a.c'7 

5-11 

j.  i'2  'Z 

5-  12a 

-  er/° 

5-  18a 

kt  ** 

6,76 

1.02 

5-11 

9/2 

5-12b 

5-  18b 

k2  * 

2  ('<7 

0  2  8 

5-11 

O.  02 

5-  13a 

-212° 

5-  19a 

k2  ** 

733 

2.66 

5-11 

c .  70. O 

5-  13b 

-7S* 

5-  19b 

k  * 

3 

7.  */ 

7.60 

5-11 

y\  t  s“" 

& .  /  »./ 

5-14 

-233* 

5-20 

■*  Relative  to  G  **  Relative  to  G  ***  0)  =  /c/  G 
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GENERAL  DATA  (for  all  pure  third  order) 


oo  (or  Q) 

M 

dM/duo 

N 

dN/dtu 

0 

1.0 

0 

0° 

— 

(JO 

M 

0 

_ 

_ 

P 

P 

00 

M 

0 

_ 

_ 

n 

m 

CO 

0 

0 

-270° 

0 

Fig  5-22  (page  69  )  compares  the  frequency  response  obtained  from 
this  rapid  approximation  method  to  the  actual  M  and  N  contours 
obtained  from  a  digital  computer  program  of  the  general  M  and  N 
equations  . 


-68- 


6^ 


70 


Ilir  r.-j.Tl  LL. 


S: 


St 


m 


t.r.  t 


jfiv- 


■St 


atw 


ipr.— 


:t± 


:# 


±H± 


m 


ROBE 


1:4  F 


"hi  : 

phi 

Hr  ‘ 

,*  |  |  1 

1  j-fif 

IT!: 

Tv 

(Tut 

tlU 

UfL 

;-Ht 

Lp  T 

::W\ 

,Tret. 

- 

.siiiti 


m 

I 


I  mm 


-44 


11 


:.! 


ss 


ts 


::T- 


m 


s 


m 


ttti 


ffifi 


Is 


t# 


$8 


It 


B- 


II 


i 


rt — !  s- 


ST 


I3± 


i 


»■ 


m 


KB 

s 


E±' 


m 


* 


ff 


r  4: 


I. 


mm 


5 


±±.  * 


$ 


IB 


t:S 


IB 


T~~ 


SS 


Ita 


m 


Btl 


:$ 


I 


ftq ' 


8S- 


s 


tst 


c±z: 


It' 


W'- 


a 


L-l 

I 


I 


ftt 


m 


rv  'i  xi 


Silt: 


| 

3 


£5 


ti 


1: 


Pig.  5-3>a 

for  Pure  Third  Order  System 

I  ;  1  Ijj  T?  flT  i  ;  l\  \  i 


_ |i.t 


m 

tit 


S 


44 


i 


ii 


m 


tr 


i-* 


■| 


m 


....  f  rl  ~ 

ESpS 


I 


X 


-Jxxnit 

itHW 


fiM 

purr 


it: 


tg 


_ 


:r.t p 


tat 


fill  ■fil  -fi 


ifi 


‘xtt 


m 


.^TU: 


Sfe 


It: 


13 


S3 


B 


BIB 


M 


Tffi 

fisr 


1  iS-JS 

iiss 


lit 


Lit 


■  4  ‘  OS, 

lUmmi 


if 
61 


is 


m lit 


s 


S 


SS  SH 

88 


II 


II 


5 


§E 

Tt" 

It  - 

,  8- 

tftt: 

illi. , 

Ii 

•f 

|bi 

]j*rf  ■ 

■fit; 

’  v-H 

.  hit : 

Tl, 

*  'ii* 

T 

"  1th- ~ 

# 

iitufi 

:l£l. 

- 1  ryf  »■  ; 

--TIT 

-1-1 


jfi 


mm. 


i]6 


y^:\ 


1 


m 


ii 


it 


mt 


st 


1.ii 


it  t 


m 


38 


-n 


iStlffl 


m 


Hffi+s 


ttX 


:s 


1 


Bti 


S 


m 


Tv 


PS 


St 


13833 


13 


~£z 


m 


Si 


<5 


Ipfi 


.SS 


It 


I! 


|SH 

ifTilLriiS.il 


It 


ffiXff! 


m 


msismst 


aiiWim 

iH1* 


tiSrt 


lit 


lit; 


4± 


T± 


m 


:tit! 


tti, 

Pit 


fw 


m 


SH 


-T-i  j 


ffi; 


tr 


I44.I 


rl 

t  xf 

TT 

T 

a 


tir 


1 


-4 

ffli 


; 


3a 


i 


i 

i 


3S 


8 


X 


1 

|±t 


il 


88'  3£ 

llL 

1 

Sat 

Us 


8S 


Fig.  5-6a 

Mp  for  Pure  Third  Order  System 
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6  CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  were  made  based  upon  the  investigation 
of  the  three  classes  of  servomechanisms  studied  in  this  thesis: 

a,  All  pure  second  order  servomechanisms  with  the  same  relative 
complex  pole  location  will  have  the  same  frequency  response. 

b  For  the  pure  second  order  servomechanism,  the  parameter  loci 
plotted  are  made  applicable  to  any  location  of  complex  poles  in  the 
left  half  of  the  s~plane  by  changing  the  scales  on  the  complex  plane 
and  on  the  parameter  frequency  loci. 

c.  The  frequency  response  for  servomechanisms  with  three 
singularities  can  be  expressed  as  a  function  of  the  referred  real  and 
imaginary  parts  of  the  complex  poles  and  the  referred  frequency  where: 

1  The  referred  real  part  of  the  complex  poles  is  the 
actual  real  part  of  the  complex  poles  divided  by  the 
magnitude  of  the  real  pole  or  zero. 

2.  The  referred  imaginary  part  of  the  complex  poles  is  the 
positive  imaginary  part  of  the  actual  complex  pole 
divided  by  the  magnitude  of  the  real  pole  or  zero. 

3»  The  referred  frequency  is  the  actual  frequency  divided 
by  the  magnitude  of  the  real  pole  or  zero. 

d.  All  servomechanisms  with  three  singularities  which  have  the 
same  re lat i ve  pol e -zero  locations  will  have  the  same  frequency  response 
relative  to  the  referred  frequency. 

e.  Curves  for  determining  values  of  several  optional  "fill-in" 

parameters  were  included  in  addition  to  the  normally  accepted  parameters 

of  M  ,  tt)  and  0),  . 

p(A)  p  bw 

f.  An  accurate  graphical  approximation  of  the  frequency  response 
of  a  servomechanism  can  be  rapidly  obtained  using  the  curves  produced 
in  this  thesis. 


1 t  is  recommended  that  the  following  areas  be  considered  for  further 
investigation  using  the  technique  presented  in  this  paper; 

a.  Servomechanisms  with  more  than  three  singularities,  and 

b.  Systems  other  than  servomechanisms. 
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8 .  Appendix  A 


The  Control  Data  Corporation  Model  1604  digital  computer 
important  characteristics  are  summarized  as  follows: 
h  All  transistorized 
2.  48  bit  word  length 

3  Stored  program 

4  Single-address  logic,  two  instructions  per  48-bit  word 
5.  Indirect  addressing  feature 

6  32,768  48-bit  words  of  magnetic  core  storage: 

(a)  Stored  in  two  independent  16,384  word  banks, 
alternate  phased 

(b)  4  8  micro  -  seconds  effective  cycle  time 
(representative  program). 

(c)  6.4  microseconds  total  cycle  time. 
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